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Abstract 

Background: The p-lactamase {bio) locus, which confers resistance to penicillins only, may control the transcription of 
mecA, the central element of methicillin resistance, which is embedded in a polymorphic heterelogous chromosomal 
cassette (the SCCmec element). In order to assess the eventual correlation between bio allotypes and genetic lineages, 
SCCmec types and/or p-lactam resistance phenotypes, the allelic variation on the bio locus was evaluated in a 
representative collection of 54 international epidemic methicillin-resistant Stophylococcus aureus (MRSA) clinical strains 
and, for comparative purposes, also in 24 diverse methicillin-susceptible 5. aureus (MSSA) strains. 

Results: Internal fragments of bloZ (the p-lactamase structural gene) were sequenced for all strains. A subset of 
strains, representative of blaZ allotypes, was further characterized by sequencing of internal fragments of the blaZ 
transcriptional regulators, blol and blaRl. Thirteen allotypes for bloZ, nine for blol and 12 for blaRl were found. In a 
total of 121 unique single-nucleotide polymorphisms (SNP) detected, no frameshift mutations were identified and 
only one nonsense mutation within blaZ was found in a MRSA strain. On average, blaZ alleles were more 
polymorphic among MSSA than in MRSA (14.7 vs 1 1.4 SNP/allele). Overall, blaRl was the most polymorphic gene 
with an average of 24.8 SNP/allele. No correlation could be established between bio allotypes and genetic lineages, 
SCCmec types and/or p-lactam resistance phenotypes. In order to estimate the selection pressure acting on the bio 
locus, the average dN/dS values were computed. In the three genes and in both collections dN/dS ratios were 
significantly below 1. 

Conclusions: The data strongly suggests the existence of a purifying selection to maintain the bio locus fully 
functional even on MRSA strains. Although, this is in agreement with the notion that in most clinical MRSA strains 
mecA gene is under the control of the bio regulatory genes, these findings also suggest that the apparently 
redundant function of blaZ gene for the MRSA resistant phenotype is still important for these strains. In addition, 
the data shows that the sensor-inducer blaRl is the primary target for the accumulation of mutations in the bio 
locus, presumably to modulate the response to the presence of p-lactam antibiotic. 

Keywords: |3-lactamase p-lactam resistance, allelic variation, MSSA, MRSA, mecA stabilization 



Background 

Staphylococcus aureus is a leading cause of nosocomial 
infections and has recently emerged as a community 
acquired pathogen [1-3]. S. aureus is also a paradigm of 
adaptive power to antimicrobial chemotherapy, able to 
develop resistance to virtually all classes of antibiotics 
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[4], The acquisition of resistance to P-lactam antibiotics 
is particularly relevant in clinical terms. Although p -lac- 
tams (i.e. penicillin G) were the first class of large-spec- 
trum antibiotics to be introduced into clinical practice, 
they are still the most widely used due to their high 
effectiveness, low cost, ease of delivery and minimal side 
effects [5]. 

In response to P-lactam chemotherapy, S. aureus has 
sequentially acquired two resistance genes: first blaZ, 
which codes for a p-lactamase and confers resistance to 
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penicillins only, and then rnecA, which codes for an extra 
penicillin-binding protein (PBP2a) with reduced affinity 
for virtually all P -lactams [6,7]. The transcription of both 
resistance genes may be controlled by homologous two- 
component systems consisting on a sensor-inducer 
(BlaRl and MecRl) and a repressor (Blal and Mecl). 
Interestingly, in spite of the cross-resistance to virtually 
all p -lactams provided by rnecA, the great majority (> 
95%) of contemporary MRS A are still positive for the P- 
lactamase locus [8]. Moreover, the regulators of blaZ, 
BlaRl and Blal, can efficiently induce rnecA transcription 
and, do it faster than the "natural" rnecA regulators, 
MecRl and Mecl [9,10]. In addition, since many MRS A 
strains do not have functional mecI-mecRl genes due to 
polymorphisms in the rnecA regulatory region [11], the 
rnecA transcription is presumably under the control of 
the blal-blaRl genes only. In line with these observa- 
tions, the presence of the blaZ locus has been shown to 
promote rnecA acquisition and stabilization [12,13]. 

In S. aureus, the p -lactamase genes may be located in 
a plasmid or mobilized into the chromosome by trans- 
poson Tn552 [14]. In contrast to the diversity of P -lacta- 
mase genes found in gram-negative bacteria, all 
staphylococcal enzymes studied so far are molecular 
class A serine p -lactamases placed in functional group 
2a [8]. The mature form of the enzyme has a molecular 
mass of 30 kDa, contains 257 amino acids, and is 
secreted extracellularly [15]. In 1965, Richmond pro- 
posed the subdivision of staphylococcal P -lactamases in 
four serotypes [16], but the structural basis of the dis- 
tinction between types is still uncertain and no clear 
relationship between sequence and serotype was found 
[17]. Interestingly, serotypes were shown to have specific 
geographic distributions [8], which may suggest a rela- 
tionship between bla-type and genetic lineage. Recently, 
Olsen et al have studied the allelic variation of the blaZ 
gene among several staphylococcal species and 11 BlaZ 
protein types were identified [14]. The multiple- 
sequence alignment of those sequence types suggest a 
separate evolution for plasmid- and chromosomally- 
encoded blaZ and a very low frequency for exchange of 
the P -lactamase locus between strains and species. 

In evolutionary terms, MRSA may be regarded as a 
recent sub-branch of the S. aureus population which has 
acquired the heterelogous chromosomal cassette con- 
taining the rnecA gene - the SCCmec element [18]. 
Molecular epidemiology studies on large collections of 
MRSA isolates have clearly shown that MRSA has a 
strong clonal structure and that very few lineages, 
defined by specific macro-restriction patterns of chro- 
mosomal DNA and/or multi-locus sequence types, 
account for the great proportion of MRSA infections 
worldwide [19,20]. The clonal structure of MRSA popu- 
lation may result from a "host barrier" for the rnecA 



acquisition, which restricts the number of acquisitions 
to few more permissive lineages [13,21] and/or from the 
clonal expansion of previously highly epidemic (MSSA) 
lineages, which have acquired the rnecA gene. Recent 
data based on comparative genomics of MRSA lineages 
[22-24] supports both mechanisms as it seems that, 
within the same genetic (epidemic) lineage, SCCmec 
acquisitions may occur continuously at the local level. 

In spite of the several lines of evidence suggesting an 
important role of the bla locus in the acquisition, stabili- 
zation and regulation of the rnecA gene, the variability of 
bla genes at the sequence level has never been evaluated 
among pandemic MRSA lineages. The present study was 
conducted in order to evaluate the allelic variability of P- 
lactamase locus in a representative collection of interna- 
tionally epidemic MRSA clones and also, for comparative 
purposes, in a diverse collection of methicillin-susceptible 
S. aureus strains (MSSA), in an attempt to make evolu- 
tionary correlations between p -lactamase allotypes and 
P -lactam resistance phenotypes (i.e. MRSA vs MSSA), 
SCCmec types and/or genetic lineages. 

Methods 

Strain collection 

S. aureus strains used in the present study are listed in 
Tables 1 (MRSA) and 2 (MSSA). All strains have been 
previously assigned to genetic lineages by Pulse-field gel 
electrophoresis (PFGE), multi-locus sequence typing 
(MLST) and protein A sequence typing (spa typing). 
MRSA strains have been additionally characterized in 
terms of their SCCmec types. The presence of a func- 
tional p -lactamase locus was confirmed by nitrocefin 
disks (Sigma) for all strains, in the presence and absence 
of an inducer (oxacillin at 0.05 mg/L). 

Media and growth conditions 

Strains were grown overnight at 37°C on tryptic soy agar 
or tryptic soy broth under aerobic conditions. 

DNA isolation 

Total DNA was prepared using the Wizard genomic 
DNA preparation kit (Promega, Madison, WI, USA), 
according to the manufacturer's recommendations, 
except for the addition of lysostaphin at 0.5 mg/mL and 
RNase at 0.3 mg/mL for the lysis step. 

DNA amplification and sequencing 

The allelic variation on the p -lactamase locus was evalu- 
ated by sequencing internal fragments of blaZ and its 
transcriptional regulators, blal and blaRl, amplified by 
PCR. Based on the available sequence at GenBank 
(accession number: X52734) for Tn552 of S. aureus, 
three pairs of primers were designed as follows (5' — > 
3'): blaZ Fl, GAT AAG AGA TTT GCC TAT GC; blaZ 
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Table 1 Characteristics and bla locus allotypes of the MRSA strains used in this study 



Clonal complex a) 


MLST (ST) 


SCCmec type 


Strain 


Isolation origin 


Isolation year 


bla locus alleles 


Ref. 














blaZ 


blal 


blaRl 






247 


I 


E2125 


Denmark 


1964 




ND 


ND 


[30,31] 




247 


IA 


HPV107 


Portugal 


1992 




ND 


ND 


[30,32] 




247 


IA 


BK1953 


USA 


1995 




Kin 

IN u 


Kin 

IN u 


[30,33] 




250 


| 


COL 


UK 


1965 








[30] 




250 


| 


BK793 


Eavot 


1961 




i 
i 


i 
i 


[30,34] 




250 


IA 


PER34 


Spain 


1989 




Kin 


Kin 

INU 


[30,35] 




239 


III 


ANS46 


Australia 


1982 




1 
1 


1 
1 


[30,36] 




239 


IIIA 


HU25 


Brazil 


1993 




Kin 

IN U 


Kin 

IN U 


[30,37] 




239 


IIIA 


HUSA304 


Hungary 


1993 




Kin 

IN U 


Kin 

IN U 


[30,38] 


8 


239 


IIIA 


BK2421 


USA 


1996 




ND 


ND 


[30,34] 




8 


IVa 


USA300 


USA 


1 995-2003 




4 


9 


[39,40] 




8 


IVa 


USA500 


USA 


1995-2003 


1 


ND 


ND 


[39,40] 




8 


IVd 


BK2529 


USA 


1996 


1 


ND 


ND 


[30,34,39] 




8 


IVd 


Bargll17 


USA 


1996 


1 


ND 


ND 


[30,39,41] 




8 


IVE 


AR43/3330.1 


Ireland 


1988-2002 


] 


1 


1 


[42] 




8 


IVh 


GRE120 


Greece 


1993 


1 


ND 


ND 


[39,43] 




72 


IVa 


USA700 


USA 


1995-2003 






4 


[39,40] 




254 


IVh 


HAR36 


Greece 


2002 


1 


1 


1 


[39,44] 




770 


IVb 


8/6-3P 


Chicago 


1996 




3 


6 


[45] 




5 


I 


HAR21 


Finland 


2002 




ND 


ND 


[44,46] 




5 


l-VAR 


PL72 


Poland 


1991 




Kin 

IN U 


Kin 

IN U 


[30,47] 




5 


|| 


N315 


Japan 


1982 


Q 
O 


A 


Q 


[18] 




5 


|| 


JP1 


Japan 


1987 


8 


4 


9 


[30,48] 




5 


II 


BK2464 


USA 


1996 


4 


U 


2 


[30,49] 




5 


|| 


USA100 


USA 


1 995-2003 


J 


J 


o 


[40,46] 




5 


IVa 


BM18 


USA 


1989 


A 


O 


Z 


[30,39,50] 




5 


IVa 


FFP31 1 


Portugal 


1996 


1 1 


] 


7 


[39,51] 


5 


5 


IVa 


HSA49 


Portugal 


1993 


1 1 


ND 


ND 


[39,51] 




5 


IVa 


HSA74 


Portugal 


1993 


5 


3 


3 


[39,51] 




5 


IVc 


Q2314 


Dallas 


1996 


3 


ND 


ND 


[52] 




5 


IVc 


USA800 


USA 


1995-2003 


1 


ND 


ND 


[39,40] 




5 


IVc 


ARG9 


Argentina 


1996 


1 1 


7 


7 


[39,51] 




5 


IVd 


JCSC4469 


Japan 


1982 


1 


1 


1 


[53] 




5 


IVg 


M03-68 


Korea 


2003 


3 


ND 


ND 


[54] 




5 


IV NT 


COB3 


Colombia 


1996 


6 


5 


10 


[30,39,55] 




5 


VI 


HDE288 


Portugal 


1996 


10 


_ 


5 


[51,56] 




5 


VI 


HUC136 


Portugal 


1995 


10 


_ 


5 


[51,56] 




228 


1 


HAR40 


Belgium 


1995 


2 


3 


1 


[44,46] 




30 


IVc 


DEN2946 


Denmark 


2001 


1 


1 


1 


[39,57] 


30 


30 


IVc 


DEN2294 


Denmark 


2001 


1 


ND 


ND 


[39,57] 




36 


II 


USA200 


USA 


1995-2003 


1 


1 


1 


[40,46] 




36 


II 


HAR24 


Finland 


2002 


1 


ND 


ND 


[46,47] 




22 


IVh 


HAR22 


Finland 


2002 


9 




4 


[39,44] 


22 


22 


IVh 


HGSA146 


Portugal 


2003 


9 


ND 


ND 


[39,58] 




22 


IVh 


HGSA163 


Portugal 


2003 


9 


ND 


ND 


[39,58] 




45 


II 


USA600 


USA 


1995-2003 


7 


4 


9 


[40,46] 


45 


45 


IVa 


HAR38 


Belgium 


1995 


6 


2 


8 


[39,44] 




45 


V 


WIS 


Australia 


1995 


8 


ND 


ND 


[59] 




256 


IVa 


CA05 


Chicago 


1999 


8 


4 


9 


[45] 
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Table 1 Characteristics and bla locus allotypes of the MRSA strains used in this study (Continued) 



1 


1 


IVa 


MW2 


USA 


1998 


6 


2 


10 


[60] 




1 


IVa 


USA400 


USA 


1995-2003 


6 


2 


10 


[39,40] 


80 


80 


IVc 


DEN2949 


Denmark 


2001 


5 


3 


3 


[39,57] 




80 


IVc 


DEN114 


Denmark 


2001 


5 


1 


3 


[39,57] 


Singleton 


377 


V 


HT0184 


Greece 


2005 


6 


2 


10 


[61] 




377 


V 


HT0826 


France 


2003 


6 


ND 


ND 


[61] 



a) Clonal complexes were determined using the E-burst software http://saureus.mlst.net/, last accessed on June 04, 2009. 
ND, not determined, -, negative PCR amplification. 



Rl, GCA TAT GTT ATT GCT TGA CC; blal Fl, GCA 
AGT TGA AAT ATC TAT GG; blal Rl, GAA AGG 
ATC CAT TTT CTG TAC ACT CTC ATC; blaRl Fl, 
CAT GAC AAT GAA GTA GAA GC; and blaRl Rl, 
CTT ATG ATT CCA TGA CAT ACG. The predicted 
amplicon sizes were 533 bp for blaZ, 484 bp for blal 
and 537 bp for blaRl, PCR was performed in a Tl 
Thermocicler (Biometra) with the following conditions: 
94°C for 4 min; 30 cycles of 94°C for 30 s, 55°C for 30 s 
and 72°C for 1 min; and a final extension at 72°C for 10 



min. In each reaction (final volume of 50 uL), 5 ng of 
chromosomal template, 2.5 U of GoTaq flexi DNA poly- 
merase (Promega), lx Colorless GoTaq flexi buffer (Pro- 
mega), 2.5 mM MgCl2 (Promega), 40 uM of each 
deoxynucleoside triphosphate (dNTPs mixture, Bioron) 
and 20 pmol of the forward and reverse primers were 
used. The amplified fragments were purified using a mix 
of Exonuclease and SAP enzymes. Sequencing of both 
strands was performed by Macrogen http://www.macro- 
gen.com or STAB Vida http://www.stabvida.com. 



Table 2 Characteristics and bla locus allotypes of MSSA strains used in this study a) 


Clonal complex b) 


MLST (ST) 


PFGE type 


Strain 


Origin 


Isolation date 




bla locus alleles 














blaZ 


blal 


blaRl 




1 


G 


IPOP38 


Portuga 


2001 


6 


2 


10 


1 


188 


L 


IPOP58 


Portuga 


2001 


6 


2 


10 




573 


M 


HSJ109 


Portuga 


1995 


6 


2 


10 


5 


5 


C 


HSA29 


Portuga 


1992-1993 


11 


4 


7 




5 


C 


IPOP41 


Portuga 


2001 


6 


3 


6 


8 


8 


J 


IPOP65 


Portuga 


2001 


8 


2 


ND 




615 


E 


IPOP32 


Portuga 


2001 


9 


1 


4 


9 


9 


D 


HSJ122 


Portuga 


1995 


12 


1 


12 


10 


10 


Q 


DCC300 


Portuga 


1996-1997 


9 


1 


5 


12 


12 


X 


HSJ130 


Portuga 


1995 


3 


3 


6 




12 


X 


Draftees728 


Portuga 


1996-1997 


1 


1 


1 


15 


15 


K 


HSA9 


Portuga 


1992-1993 


6 


9 


ND 


20 


20 


N 


HSA47 


Portuga 


1992-1993 


6 


8 


11 


22 


22 


T 


Draftees721 


Portuga 


1996-1997 


6 


3 


5 


25 


25 


S 


HSA76 


Portuga 


1992-1993 1 1 1 




30 


A 


IPOP37 


Portuga 


2001 


13 


1 


1 


30 


34 


B 


IPOP24 


Portuga 


2001 


6 


ND 


ND 




34 


B 


IPOP34 


Portuga 


2001 


1 


1 


ND 




NA 


B 


IPOP26 


Portuga 


2001 


1 


ND 


ND 


45 


45 


H 


HSA19 


Portuga 


1992-1993 


6 


2 


10 




45 


H 


IPOP56 


Portuga 


2001 


6 


ND 


ND 


97 


97 


P 


IPOP50 


Portuga 


2001 


6 


ND 


ND 


121 


121 


F 


IPOP44 


Portuga 


2001 


10 


1 


5 


Singleton 


580 


R 


DCC1185 


Portuga 


1996-1997 1 1 1 



a) MSSA strains have been previously characterized by PFGE and MLST [62]. 

b) Clonal complexes were determined using the E-burst software http://saureus.mlst.net/eburst/database.asp, last accessed on June 04, 2009. 
NA, not available; ND, not determined. 
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DNA sequences analysis and phylogenetic tree 
reconstruction 

DNA sequencing raw data analysis and multi-sequence 
alignments were performed using the DNA Star soft- 
ware package (Lasergene). For the multi-sequence align- 
ments, the Clustal W algorithm was used. In order to 
maximize sequence reads, raw sequences for blaZ and 
blaRl were trimmed immediately after the primer 
sequences keeping the reading frame. As the reverse pri- 
mer for blal (Blal Rl) is located outside of the coding 
region, the 3' end of the sequence was trimmed at the 
end of the coding region. For each gene, allotypes were 
defined taking as reference the extant sequences of the 
bla locus of Tn552, which were assigned to allotype 1. 

Phylogenetic and molecular evolutionary analyses were 
conducted using MEGA version 4 [25] and the resultant 
phylogenetic trees were obtained using the neighbour- 
joining (NJ) method with bootstrap analysis using 1000 
replicates. In order to evaluate the diversity of the bla 
locus, the Simpsons indexes of diversity (SID) were cal- 
culated [26,27] for each locus using the online tool 
available at http://www.comparingpartitions.info. To 
estimate selection pressure acting on the bla locus, we 
computed the dN/dS ratios for the three genes. The 
dN/dS ratios were computed for all pairs of alleles with 
more than 1% substitutions, in order to give an estimate 
of the divergence of the alleles while excluding those 
pairs that, being too similar, would give anomalous dN/ 
dS ratios. The dN/dS ratios were computed by Model 
Averaging, as described in [28] and implemented in the 
KaKs_Calculator application [29]. This approach fits a 
set of models by maximum likelihood and then com- 
putes the weighted average of the models using a sec- 
ond-order Akaike Information Criterion (AICC). 

Nucleotide sequence accession numbers 

All nucleotide sequences determined in this study were 
deposited in Genbank under accession numbers 
GQ980053-GQ980139 (blaZ alleles), GQ980140- 
GQ980187 (blal alleles) and GQ980188-GQ980236 
(blaRl alleles). 

Results 

The allelic variation in the p -lactamase locus (bla) was 
evaluated by sequencing internal fragments of blaZ, blal 
and blaRl genes in a representative collection of inter- 
national epidemic MRSA clones and also, for compara- 
tive purposes, in a diverse collection of MSSA strains. 

blaZ allelic variability 

Thirteen different blaZ allotypes were identified within 
our collection, which comprised 54 MRSA and 24 
MSSA (Tables 1 and 2, respectively). Although seven 
alleles were common to MRSA and MSSA strains, we 



found four alleles present in MRSA strains only and two 
present in MSSA strains only. Moreover, the relative fre- 
quencies of each allele were different among MRSA and 
MSSA strains (Table 3); for instance, blaZ allotype 1 
was dominant in MRSA strains accounting for 43% (23 
out of 54) of the isolates whereas in MSSA it accounted 
for 21% (5 out of 24) of the isolates, and blaZ allotype 6 
was present in 11% (6 out 54) of MRSA but was domi- 
nant among MSSA accounting for 46% (11 out 24) of 
the isolates. The diversity of blaZ gene as measured by 
the Simpson index of diversity (SID) was higher for the 
MRSA collection than for MSSA, although not statisti- 
cally significant due to the partial overlapping of the 
confidence intervals (SID = 79.18, 95%CI 69.6-88.8 vs 
SID = 76.09, 95%CI 61.3-90.9, respectively) - see Table 
4. Within the length of blaZ region analyzed (492 
nucleotides), we detected 43 unique single-nucleotide 
polymorphisms (SNP) and on average, each blaZ allele 
has 12.4 SNP comparing to the prototype blaZ sequence 
of Tn552 (allele 1) - see Tables 3 and 4. Overall, blaZ 
alleles were more variable in MSSA than in MRSA (14.7 
and 11.4 SNP/allele, respectively). As illustrated by the 
allelic frequency distribution per MRSA lineage (Figure 
1) or the cluster tree of the thirteen blaZ alleles found 
in our collections (Figure 2), there is no clustering 
according to genetic lineages, as defined by MLST 
sequence type and SCCmec type, or MSSA/MRSA phe- 
notype; i.e. the same allele could be detected in different 
genetic lineages or among MRSA and MSSA, and the 
same lineage could be characterized by several alleles. In 
addition, there was also no clear clustering of blaZ allo- 
types according to geographic origin or isolation date of 
the MRSA isolates (see Table 1). 

The BlaZ variability in the MRSA and MSSA strains at 
the protein level was evaluated by comparison of the 
deduced amino acid sequence of all alleles against the 
deduced amino acid sequence for the BlaZ of Tn552. 
Overall, the deduced amino acid sequences of blaZ 
alleles from the MRSA and MSSA strains revealed on 
average 5.8 silent mutations, 1.8 conservative missense 
mutations and 4 non-conservative missense mutations 
per allotype (see Tables 3 and 4). For MRSA strain 
HAR40, a nonsense mutation at Gln76 was detected 
which presumably originates a non-functional truncated 
BlaZ protein. As this strain was positive for the nitroce- 
fin test, the DNA extraction and the blaZ sequencing 
were repeated and the nonsense mutation was con- 
firmed. No frameshift mutations were found in blaZ 
allotypes. 

Allelic variability of blaZ regulatory genes 

Based on the blaZ variability analysis, we selected 51 
representative strains to further characterize the variabil- 
ity in the blaZ regulatory genes, blal and blaRl, Some 
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Table 3 Characteristics of bla locus alleles 



Gene 


Allele No. 


Frequency 


SNP C) 




Amino acid substitutions 








MRSA a) 


MSSA b) 




Silent 


Conservative 


Missense 


Nonsense 




1 


U.4o 


n 1 1 
U.z I 


0 


0 


U 


U 


0 




2 


U.Uz 


U 


1 


o 


n 

U 


i 
l 


1 




3 


U.U/ 


U.U4 


9 


4 


Z 


z 


o 




4 


0 04 


o 


9 


4 


2 


3 


o 




5 


0.06 


o 


7 


2 


2 


3 


0 




6 


0.1 1 


0.46 


13 


8 


2 


3 


0 


blaZ 


7 


0.02 


o 


12 


6 


2 


4 


0 




8 


0.10 


0.04 


11 


6 


2 


3 


0 




9 


0.07 


0.08 


20 


9 


2 


7 


0 




10 


0.04 


0.04 


19 
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a) The total number of MRSA strains whose blaZ, blal and blaRl genes were analyzed is 54, 27 and 31, respectively. 

b) The total number of MSSA strains whose blaZ, blal and blaRl genes were analyzed is 24, 20 and 17, respectively. 

c) For each allele, the SNP were counted taking as reference Tn552 bla sequences (allele 1). 



of these strains failed in the amplification of one of the 
blaZ regulatory genes (see Tables 1 and 2). 

Within the length of blal region analyzed (351 nucleo- 
tides), we detected 13 unique SNP, which account for 
the nine blal allotypes detected (see Tables 3 and 4). 
Four of the nine blal allotypes were present in both 
MRSA and MSSA, while three blal allotypes were found 
in MRSA strains only and two in MSSA only. The SID 
was higher for MRSA than for MSSA although not sta- 
tistically significant (SID = 82.1, 95%CI 74.6-89.5 vs SID 
= 74.2, 95%CI 60.5-87.9, respectively) (Table 4). On 
average, each blal allele has 3.4 SNP comparing to the 



prototype blal sequence of Tn552 (allele 1), and blal 
alleles were on average more polymorphic for MRSA 
than for MSSA (3.9 vs 2.5 SNP per allele, respectively) - 
see Tables 3 and 4. 

Within the length of blaRl region analyzed (498 
nucleotides), we detected 65 unique SNP, which account 
for the 12 blaRl allotypes detected (see Tables 3 and 4). 
Six of the 12 blaRl allotypes were present in both 
MRSA and MSSA, while four blaRl allotypes were 
unique for MRSA strains and two were characteristic of 
MSSA strains. The SID values were virtually identical 
for both MRSA and MSSA (SID = 88.8, 95%CI 83.2- 
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Table 4 Comparative analysis of the allelic variation in bla locus for MRS A and MSSA strains 

No. No. Simpson's Unique SNP/allele Mutations per allele (average) dN/dS 

isolates alleles SNP (average) 



ID CI Silent Conserv. Missense Nonsense Average St. 

(95%) dev. 





blaZ 


54 


1 1 


79.2 


69.6- 
88.8 


41 


1 1.4 
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1 .7 
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MSSA 
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0.08 


0.03 




blaRl 
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81.2- 
95.3 


61 
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7.8 
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0.24 


0.10 




blaZ 


78 


13 


81.1 


75.0- 
87.3 


43 


12.4 


5.8 


1.8 


4.0 


0.1 


0.20 


0.10 


All 


blal 


47 


9 


78.4 


71.0- 
85.9 


13 


3.4 


2.3 


0.1 


1.0 


0 


0.10 


0.04 




blaRl 


48 


12 


88.5 


84.0- 
93.0 


65 


24.8 


10.2 


5.3 


8.1 


0 


0.25 


0.10 



ID, index of diversity; CI, confidence interval; SNP, single-nucleotide polymorphism; Conserv., conservative; St. dev., standard deviation 



94.4 vs SID = 88.2, 95%CI 81.2-95.3, respectively) (Table 
4). On average, each blaRl allele has 24.8 SNP compar- 
ing to the prototype blaRl sequence of Tn552 (allele 1), 
with no significant differences between MRSA and 
MSSA (24.4 and 24.6 SNP/allele, respectively) - see 
Tables 3 and 4. 

In agreement with what was observed for the blaZ 
gene, the cluster trees of blal and blaRl alleles found in 
our collections also showed no clustering according to 
MSSA/MRSA phenotype or genetic lineages (Figures 3 
and 4). For those strains in which the alleles of the three 
genes were determined, we constructed a cluster tree 
with the concatenated sequences - see Figure 5. In spite 



of the relatively low number of allelic profiles, there was 
still no clear clustering of bla allotypes according to 
MSSA/MRSA phenotype or lineage, as the same allelic 
profile was present in different genetic lineages (e.g. pro- 
file 8/4/9 present in clonal complexes 5, 8 and 45) and, 
the same genetic lineage was characterized by profiles 
from different brunches (e.g. clonal cluster 8 character- 
ized by profiles 8/4/9, 1/1/1, 3/3/6, etc.). 

The Blal and BlaRl variabilities at the protein level 
in the MRSA and MSSA strains were evaluated by 
comparison of the deduced amino acid sequence of 
all alleles against the corresponding deduced amino 
acid sequences of Tn552 (see Tables 3 and 4). 
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Figure 1 blaZ allotype frequency per MRSA lineage as defined 
by MLST clonal cluster. 
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Figure 2 Cluster tree of blaZ gene allotypes found in the MRSA 
and MSSA collections. The tree was constructed with the 
neighbor-joining (NJ) method. In each branch is shown the 
corresponding bootstrap NJ values, taken over 1000 replicates, 
which assigns confidence values for the groupings in the tree. For 
each allele, it is indicated the collection(s) (MRSA or MSSA) and 
genetic lineage (clonal cluster) where it was found. 
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Figure 3 Cluster tree of blal gene allotypes found in the MRSA 
and MSSA collections. See Figure 2 legend for details. 



Overall, the deduced amino acid sequences of the 
blal alleles revealed on average 2.3 silent mutations, 
0.1 conservative missense mutations and 1 non-con- 
servative missense mutation per allotype. The 
deduced amino acid sequences of the blaRl alleles 
showed on average 10.2 silent mutations, 5.3 conser- 
vative missense mutations and 8.1 non-conservative 
missense mutations per allotype. None of the SNP 
detected within the blal or blaRl resulted in non- 
sense or frameshift mutations. 

Selection pressure acting on the bla locus 

Based on the allelic data obtained, we computed the dN/ 
dS ratios as estimates for the selective pressure acting 
on the bla locus. The dN/dS ratios were computed for 
all pairs of alleles differing more than 1%, in order to 
give an estimate of the allelic divergence, excluding the 
anomalous dN/dS ratios of those pairs being very simi- 
lar. The average of the obtained dN/dS values and 
respective standard deviations are summarized in Table 
4. The dN/dS values for the three genes in the MRSA, 
MSSA and MRSA/MSSA partitions were well below 1 
(between 0.08 and 0.25 with standard deviations 
between 0.05 and 0.1), which suggests a negative or pur- 
ifying selection acting on the bla locus. In agreement 
with the average number of SNP per allele, the dN/dS 
ratios were significantly higher for the blaRl gene (0.24 
- 0.25) and lower for meel (0.08 - 0.11). 
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Figure 4 Cluster tree of blaRl allotypes found in the MRSA and 
MSSA collections. See Figure 2 legend for details. 
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Figure 5 Cluster tree of the concatenated blaZ-blaRl -blal 
sequences found in the MRSA and MSSA collections. See Figure 
2 legend for details. 



Discussion 

The rationale for this study comes from several observa- 
tions strongly suggesting a role of bla genes in the 
acquisition, stabilization and regulation of mecA gene, 
the central element of "broad-spectrum" p -lactam resis- 
tance characteristic of MRSA strains. The purpose of 
this study was to evaluate the allelic variability of the 
bla locus in a representative collection of international 
epidemic MRSA clones and also, for comparative pur- 
poses, in a diverse collection of MSSA strains, in an 
attempt to establish evolutionary correlations between 
bla allotypes and p -lactam resistance phenotypes (i.e. 
between MRSA and MSSA), SCCmec types (i.e. poly- 
morphisms in the mecA regulatory locus) and/or genetic 
lineages. 

MRSA lineages are much less diverse than MSSA 
lineages in terms of their genome content, a conse- 
quence of their more recent evolutionary history [19,20] 
and, apparently, also due to some "host barrier" to the 
SCCmec acquisition [13]. These differences in genetic 
background variability were well illustrated in our col- 
lections since the international MRSA collection com- 
prised eight lineages as defined by MLST clonal 
complexes, whereas in the smaller and local MSSA col- 
lection 15 lineages were represented. 

In contrast to the genetic background diversity, we 
could not detect significant differences between MSSA 
and MRSA in terms of the bla locus allelic variability. 
Actually, there were disparate subtle differences in 
terms of number of allotypes and number of point 
mutations per allotype: e.g. 11 vs 9 blaZ allotypes and 
11.4 vs 14.7 SNP/allele in MRSA and MSSA, respec- 
tively. These subtle differences may reflect the more 
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ancient evolutionary history of MSSA or a selective 
pressure to improve the bla locus activity in these 
strains. That is to say, although fewer bla types have 
been retained by the natural selection in MSSA, on 
average, these allotypes seem to have accumulated more 
adaptive mutations, in comparison to MRSA strains. In 
particular for blaZ, for which differences in terms of 
number of alleles and SNP/allele were more significant, 
the presence of the alternative p -lactam resistance 
mechanism mediated by the mecA gene in MRSA 
strains might have allowed a release in the selective 
pressure to keep blaZ with optimal activity, in contrast 
to MSSA, which rely only on blaZ-mediated resistance 
to p -lactams. 

No correlation could be established between bla allo- 
types and strain backgrounds, P -lactam resistance phe- 
notypes, strain origin and/or isolation dates, indicating 
that bla genes have evolved independently from S. aur- 
eus clonal lineages. This is particularly striking for 
MRSA strains, which have a very strong clonal struc- 
ture. These observations may be explained either by dif- 
ferences in evolutionary clock speeds between the 
genetic background and the bla locus or may result 
from the horizontal transfer of bla genes between differ- 
ent lineages, which are usually integrated in mobile ele- 
ments (plasmids and composite transposons). 
Interestingly, based on the characterization of a collec- 
tion of several staphylococcal species, Olsen et al, sug- 
gested that there is little exchange of bla genes between 
strains or species [14], which somehow contradicts our 
findings. In our study, the most parsimony explanation 
for the presence of the same bla type in different 
genetic lineages either MRSA or MSSA or the presence 
of several bla types in the same lineage, is indeed a high 
frequency for the horizontal transfer of bla genes across 
S. aureus clonal clusters. 

In spite of the lack of evolutionary links between bla 
allotypes and genetic lineages, our data strongly suggests 
a selective pressure to keep the bla locus fully func- 
tional, as illustrated by the calculated average dN/dS 
values well below 1. This observation is valid even on 
MRSA for which one could expect the accumulation of 
nonsense or frameshift mutations that would render the 
bla locus non-functional, due to presence of the mecA 
gene. Actually, the majority of the mutational events 
detected in this study were either silent or neutral muta- 
tions, being the blaRl the gene with the highest muta- 
tional rate and the blal the one with the lowest. The 
increased allelic variability detected for blaRl (in terms 
of number of alleles, Simpson's index of diversity, aver- 
age SNP/allele, and dN/dS values) may suggest that this 
sensor-inducer gene is the primary target for the evolu- 
tionary adaptive mechanisms in the bla locus, presum- 
ably to improve the induction efficiency of blaZ 



expression or even mecA expression, in the case of 
MRSA strains with no functional mecI-mecRl regulatory 
system. In contrast, the relatively lower variability of the 
much smaller blal gene, may suggest a fine-tuned 
repressor activity and a selective pressure to maintain 
the repressor activity; i.e to maintain the blaZ expres- 
sion inducible. 

Despite the cross-resistance to virtually all p -lactam 
antibiotics provided by mecA, most contemporary 
MRSA strains still carry, besides the SCCmec element, 
the p -lactamase locus. This might be due to the fact 
that not enough time has elapsed since the mecA acqui- 
sition for MRSA strains start loosing the bla genes, 
because there is a little or no fitness cost associated to 
the bla genes, or because these genes may be linked to 
other positively selected genes (e.g. the cadmium resis- 
tance genes present in some P-lactamase plasmids). 
Alternatively, the bla locus may be involved in the 
"domestication" of the mecA gene, as bla genes have 
been shown to stabilize the in vitro mecA acquisition 
[12,13] and efficiently control mecA transcription [9,10], 
explaining the "retention" of a functional bla regulatory 
system by most contemporary MRSA strains [8]. Inter- 
estingly, as no correlation could be established between 
bla allotypes and SCCmec types, which have polymorph- 
isms in the mecA regulatory locus, this maintenance of 
functional blal-blaRl genes seems to be independent of 
the functional status of the mecA "natural" regulators 
mecI-mecRl. 

Concerning the maintenance of a functional blaZ gene 
in MRSA strains one can speculate that, even in the pre- 
sence of mecA, it might be useful for the bacteria to 
keep blaZ as a "first-line defense" against p -lactams. In 
fact, first generation p -lactams (i.e. penicillins) are still 
widely prescribed either empirically or for the treatment 
of specific infections (e.g. streptococcal infections). 
Moreover, penicillins have also been widely used pro- 
phylactically in the livestock industry. This means that, 
both in the nosocomial and community settings, MRSA 
are still exposed to penicillins and, under these circum- 
stances, expression of P-lactamase is enough for survival 
under antibiotic pressure. From a physiological perspec- 
tive, this ability to choose between the expression of two 
resistance genes may be advantageous for the bacteria 
since the expression of P-lactamase is likely to impose a 
smaller fitness cost than the expression of PBP2a. In 
fact, besides being much smaller than PBP2a (257 vs 
668 amino acids), BlaZ is a secreted enzyme whereas 
PBP2a is a transpeptidase protein, which must be incor- 
porated into the complex cell-wall metabolism. 

Conclusion 

In this study we have evaluated the allelic variation of the 
bla locus in MRSA and MSSA clinical strains. Although 
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no correlation between bla allotypes and genetic lineages, 
SCCmec types and p -lactam resistance pheno types could 
be established, we provided evidence for the existence of 
a selective pressure to maintain the bla system fully func- 
tional even on MRSA strains and that the sensor-inducer 
gene blaRl is the primary target for the accumulation of 
adaptive mutations in the bla locus. 
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